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Q-c Abstract. Metal-dielectric multilayer metamaterials with extreme loss-anisotropy, 

o 



in which the longitudinal component of the permittivity tensor has ultra-large 
imaginary part, are proposed and designed. Diffraction-free deep subwavelength beam 
propagation and manipulation, due to the nearly flat iso- frequency contour (IFC), 
J>-^ is demonstrated in such loss-anisotropic metamaterials. It is also shown that deep 

i-G ' subwavelength beam propagation can be realized in practical multilayer structures 

with large multilayer period, when the nonlocal effect is considered. 
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1. Introduction 

Light beam with subwavelength confinement and propagation beyond diffraction limit is 
highly desirable for many optical integration applications. To enable the subwavelength 
optical confinement, optical materials supporting ultra-high wave vectors are necessary. 
Since the refractive indices of natural materials are quite limited at optical frequency, 
metamaterials with artificially engineered subwavelength meta-atoms are designed to 
exhibit ultra-high refractive indices so as to achieve large wave vectors [1, 2]. To obtain 
diffraction-free deep subwavelength beam propagation, a nearly flat IFC curve over a 
broad range in fc-space is required, so that all spatial components will propagate with 
the same phase velocity along the longitudinal direction [3]. Extremely anisotropic 
metamaterials with infinite real part of permittivity have been theoretically proposed 
to achieve flat IFC and consequently the subwavelength beam propagation without 
diffraction [4, 5]. 

Recently, metal-dielectric multilayer metamaterials with indefinite permittivity 
tensor have been utilized to demonstrate intriguing applications of negative refraction 
[6], subwavelength imaging [7], enhanced photonic density of states [8], and broadband 
light absorbers [9], together with ultra-high refractive indices for subwavelength 
optical waveguides [10] and indefinite cavities [11]. In the present work, we propose 
the concept of extreme loss-anisotropy in metal-dielectric multilayer metamaterials, 
where the longitudinal component of the permittivity tensor has ultra-large imaginary 
part. Diffraction-free deep subwavelength beam propagation and manipulation is 
demonstrated in such loss-anisotropic metamaterial, due to the nearly flat IFC. 

2. Theory and Discussion 

The inset of Fig. 1 shows the metal-dielectric multilayer structure, where titanium oxide 
(T13O5) with dielectric constant of 5.83 and silver (Ag) are chosen [12]. When the 
thickness of each layer is infinitely small, the multilayer structure can be regarded as an 
anisotropic effective medium with permittivity tensor of 

£x = fd^d + fm^m-, £y = {fd/^d + fm/^m) , (1) 

where Ed and e m are the permittivity of titanium oxide and silver, fa and f m (fd+f m = 1) 
are the filling ratios of titanium oxide and silver, respectively. The silver filling ratio 
f m is 0.45, and its permittivity is from the experimental results [13]. The dependence 
of permittivity tensor on wavelength A is shown in Fig. 1(a), which is calculated using 
Eq. (1) according to the effective medium theory. It is found that the longitudinal 
permittivity component e y shows a strong resonance at Ao = 406.1 nm (at position II), 
where lm(e y ) has a peak more than 230, while Ke(e y ) flips its sign quickly across the 
resonance from negative maximum (at position I) to positive maximum (at position III). 
Previously, ultra- large Re(e y ) (at position III) has been utilized to realize subwavelength 
beam propagation without diffraction [4, 5]. However, the behavior of ultra- high Im^) 
has not been considered before. It is intuitively thought that a large Im(e^) will 
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increase the beam propagation loss. On the contrary, it will be demonstrated that the 
extreme loss-anisotropy with ultra-large imaginary part can enable low-loss diffraction- 
free subwavelength beam propagation. For TM-polarized light with non- vanishing E x , 
E y and H z field components, the corresponding IFC is determined by 

k x / e y + k y / £ x = k , (2) 

where k x is the transverse k- vector and k y is the propagation k- vector. Figure 1(b) shows 
that an ultra-flat IFC over a large k x range is supported at the resonance wavelength 
A = 406.1 nm (position II), with e x = 1.06 + 0.098i and e y = 10.48 + 231. 70i. The 
propagation k- vector Ke(k y ) remains a constant for all different k x . While the imaginary 
part of the propagation fc-vector Im(k y ) is quite small. To understand the behavior of 
the flat IFC, Eq. (2) can be rewritten as 

/ 7 p\ £2 

K = Je x (k 2 -f\~ V^k - V^^p (3) 

where the approximation \e y \ ^> 1 has been used to derive the above formula. Since 
e x is dominated by its real part, and e y is dominated by its imaginary part at the 
resonance wavelength, Re(k y ) ~ \fe x kQ and Im(k y ) ~ ^/e^k x / (2k Im(e y )), which shows 
that Re(k y ) is independent of k x and lm(k y ) is weakly proportional to k\. Moreover, 
the propagation loss Im(k y ) is inversely proportional to lm(e y ) and the large material 
loss lm(e y ) actually enables the low-loss beam propagation. The ultra-large Im(e y ) in 
extremely loss-anisotropic metamaterial not only gives rise to an ultra-flat IFC over a 
broad k x range, but also results in ultra-small beam propagation loss. 

To further illustrate the importance of large lm(e y ), the comparisons of IFC curves 
at three different wavelengths of I, II and III are shown in Fig. 1(c), which corresponds 
to negatively maximized Ke(e y ), maximized lm(e y ) and maximized Ke(e y ), respectively. 
It indicates that the maximized lm(e y ) case results in the flattest IFC with almost zero 
curvature, while the other two IFCs show positive and negative curvatures, respectively. 
These behaviors can be clearly understood using Eq. (3), in which a purely real e y will 
contribute to Re(k y ) and influence the curvature of the IFC (a positive e y leads to a 
negative curvature and vice versa). It is clear that the maximized lm.(s y ) case with zero 
curvature will achieve diffraction-free subwavelength beam propagation. 

Next, diffraction-free deep-subwavelength beam propagation will be demonstrated 
from the numerical simulation based on finite element method (FEM). Here it is 
worthwhile to define the minimal waist size of light beam which can propagate inside 
the multilayer structure without diffraction. It is known from Eq. (2) that k 2 x je y -C k\ 
is required to obtain a nearly flat IFC, then the minimal beam waist size w m \ n turns 
out to be 2n /max(fc x ) = Aoy|&y|. Figure 2 shows the propagation of ultra-narrow 
Gaussian beams (with a waist size of 40 nm ~ O.lAo > w m in) inside the loss-anisotropic 
metamaterials with different geometries at the resonance wavelength Ao = 406.1 nm. 
The effective medium results and the realistic multilayer results are shown in Fig. 2(a- 
c) and Fig. 2(d-f), respectively. Figures 2(a) and 2(d) show that ultra-narrow Gaussian 
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beams can propagate over a long distance without any wave-front distortion. Two 
subwavelength beams with 150 nm center-to-center distance remain well-defined as the 
beams propagating across the multilayer from the bottom to the top. It is emphasized 
that the subwavelength beam confinement is entirely due to the unique loss-anisotropic 
property, and the beam path is solely determined by the launching location. This 
is distinguished from the situation in a subwavelength waveguide, where the mode is 
confined by the waveguide boundary. 

Besides the straight beam propagation, the flow of light can be flexibly modeled 
through controlling the local metamaterial properties. For instance, the beam path can 
be manipulated by gradually varying the direction of multilayers, since the direction of 
beam propagation is always vertical to the multilayer interface. The designed geometries 
for achieving 90° and 180° bending of subwavelength beams are shown in Figs. 2(b-c) 
and 2(e-f), for both the effective media and the multilayer structure. In the effective 
medium calculation, the anisotropic permittivity tensor depends on the tilted angle 
of the multilayer interface. In the local coordinate (u,v), the components of the 
permittivity tensor can still be determined by the mixing formula in Eq. (1). The 
permittivity tensor expression in the global coordinate (x,y) is related to that in the 
local coordinate as 

<_,., . I e u cos 2 9 + e v sin 2 9 (e u — e v ) sin 9 cos 9 \ . . 

I (e u — e v ) sin 9 cos 9 e u sin 2 9 + e v cos 2 9 J 

where 9 is the local tilted angle of multilayer with respect to the +x axis, e u and e v are the 
local permittivity tensor components along and normal to the multilayer, respectively. 
The simulation results indicate that the flow of light can indeed be manipulated while 
maintaining the deep-subwavelength beam confinement and diffraction-free propagation. 
For the results of 180° bending shown in Figs. 2(c) and 2(f), it is noted that the wave- 
front becomes tilted at the output section (but the energy flow is still vertical to the 
interface as a result of the flat IFC). This is due to the fact that the light traveling 
at the inner side undergoes less optical path than the light traveling at the outer side. 
The phase difference arising from the light path difference leads to the beam wave-front 
tilting. 

Figure 2 shows that the multilayer structure simulation results agree very well 
with the EMT results, indicating that the multilayer structure with period a = 20 nm 
(/m — 0.45) can represent the loss-anisotropic effective medium well. However, the 
fabrication of such thin layers is very challenging in reality (but possible [14]). It will 
be interesting to study the property of metal-dielectric multilayer with a large period 
a, where the nonlocal effect has to be taken into account [15]. The dispersion relation 
describing the realistic multilayer structure for TM-polarized light is 

cos [k y (a m + ad)] = cos(k m a m ) cos(k d a d ) - 7™ sin(fc m a m ) sm(k d a d ), (5) 

which is derived by treating the layered structure as a one-dimensional photonic crystal. 
Here 7 TM = {s d k m /e m k d + e m k d / s d k m ) / 2, k m = \ZE m k — k x , and k d = y s d k Q — k x . 
The IFCs corresponding to the realistic multilayer structure with a = 40 nm are shown 
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in Fig. 3(a). It is found that the IFC curve at Ao = 406.1 nm is no longer flat due to the 
nonlocal effect. That is to say, the effective permittivity tensor becomes strongly wave 
vector dependent, so that the permittivity components e x and e y will be functions of not 
only the frequency but also the fc-vector. The frequency corresponding to the flattest 
IFC will then be shifted. For the multilayer structure with a = 40 nm, it turns out that 
the flattest IFC occurs at Ao = 418.3 nm (which can be mathematically determined by 
finding the working wavelength with zero IFC curvature), as shown in Fig. 3(a). The 
propagation of subwavelength Gaussian beams at the two wavelengths in the multilayer 
structure with a = 40 nm are shown in Fig. 3(b). As can be expected from the IFC curves 
in Fig. 3(a), the Gaussian beams at Ao = 406.1 nm suffer strong diffraction, resulting in 
distorted beam profiles. In comparison, the beam profiles remain well-defined for the 
Gaussian beams at Ao = 418.3 nm. 

3. Conclusion 

In conclusion, extremely loss-anisotropic metamaterial is designed using metal-dielectric 
multilayer structures. The IFC corresponding to such metamaterial turns out to 
be ultra-flat over a broad fc-vector range. This unique property is then utilized 
to obtain diffraction-free deep subwavelength beam propagation. Furthermore, it is 
shown that the propagation of light beams can be manipulated flexibly by tuning the 
direction of multilayer structure. Moreover, the nonlocal effect occurring in multilayer 
structures with large multilayer period is investigated. It is found that diffraction-free 
beam propagation is still possible after taking into account the nonlocal effect. The 
current study is very attractive for many applications such as optical imaging, optical 
integration, on-chip optical communication. 
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Figure Captions 

Figure 1: (Color online) (a) The dependence of the permittivity tensor on the 
wavelength Ao for the metal-dielectric multilayer structure with metal filling ratio 
f m = 0.45. Re(e x ) has been scaled up by 10 times. The three Roman numbers (I, II, 
III) indicate the wavelengths corresponding to negatively maximized Ke(e y ), maximized 
ImfV.y) and maximized Ke(e y ), respectively. Inset shows the multilayer structure with 
TiaOs (green color) and silver (yellow color), (b) The IFC at the resonance wavelength 
Ao = 406.1 nm (at position II). ~Re(k y ) and Ira(k y ) are represented by the blue line and 
green line, respectively. Although lm(k y ) grows with transverse wave vector k x (ranges 
from 0.0476A;o to 0.2674A; as k x increases from to 10k ), it remains small due to the 
large magnitude of e y , which gives a low propagation loss. The dashed black circle is 
the IFC of air. (c) The IFCs for the three wavelengths indicated in (a). The flattest 
IFC is obtained at the resonance wavelength with maximized Im(5 y ). 

Figure 2: (Color online) Diffraction-free deep subwavelength beam propagation in 
extremely loss-anisotropic metamaterials. The distributions of magnetic field H y are 
shown in (a-c) for ideal effective medium and (d-f) for multilayer structure with period 
a = 20 nm. (a, d) Straight beam propagation, (b, e) 90° beam bending, (c, f) 180° 
beam bending. The center-to-center distance of the two beams is 150 nm in all the 
simulations. 

Figure 3: (Color online) (a) The IFCs for realistic multilayer structure with a = 
40 nm at two different wavelengths, (b) The Gaussian beams propagation at the two 
wavelengths. Diffraction-free beam propagation is achieved after taking into account 
the nonlocal effect induced wavelength shift. 
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Figure 1. Yingran He, Lei Sun, Sailing He, and Xiaodong Yang 
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Figure 2. Yingran He, Lei Sun, Sailing He, and Xiaodong Yang 
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(a) a = 40 nm (b) 400nm 
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Figure 3. Yingran He, Lei Sun, Sailing He, and Xiaodong Yang 



